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I. INTRODUCTION 

The Kerr-Taub-NUT (KTN) spacetime 0, first discovered by Demianski and Newman and later studied by 
many others (see for references and a global analysis), describes a stationary axisymmetric object with gravit- 
omagnetic monopole and dipole moments associated with nonzero values of the NUT and Kerr parameters £ and 
a respectively, and as such is a useful model for exploring gravitomagnetism. The well-known Dirac quantization 
of the magnetic monopole corresponds in the gravitomagnetic case to a condition first found by Misner 0, |^ for 
the nonrotating special case of the Taub-NUT spacetime 0] (zero Kerr parameter) that is consistent with forcing 
periodicity in the time coordinate, making the gravitomagnetic monopole Taub-NUT spacetime unphysical under 
normal circumstances by having closed timelike lines, but nonetheless interesting as a laboratory for probing other 
consequences of gravitomagnetic monopoles. As discussed by Miller 3], this same global structure is also consistent 
with the more general Kerr-Taub-NUT spacetime, which is the rotating version of this simpler spacetime, its rotation 
being associated with the nonzero gravitomagnetic dipole moment. 

The Kerr-Taub-NUT spacetime and its special cases all belong to the larger class of stationary axisymmetric 
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type D vacuum solutions of the Einstein equations found by Carter [3 for which the Hamilton-Jacobi equation for 
geodesies is separable. The stability of these spacetimes is probed by studying their perturbations by fields of various 
spin. Perturbations by massless fields have been investigated in the two limiting cases of the Kerr spacetime (pure 
gravitomagnetic dipole) ^ and the Taub-NUT spacetime (pure gravitomagnetic monopole) in a single unifying 
approach making use of the de Rham Laplacian which allows all fields of spin through 2 to be considered together. 
In this article, the same analysis is applied to the Kerr- Taub-NUT spacetime. Motivated by scattering off virtual black 
holes, Prestidge has performed this analysis for all spins except 3/2 for the vacuum C-metric, another vacuum 
type D metric representing a pair of uniformly accelerated gravitoelectric monopoles. 

This approach to perturbations builds on the pioneering work of Teukolsky 0, j done in the context of the 
Newman-Penrose formalism 0, , which partially received its mathematical foundation from Stewart and Walker 
[T5I I and was subject to some extensions by other authors [T^llTj . Teukolsky found a separable master equation whose 
eigenfunction solutions essentially solve the problem of the massless perturbations of any spin for the Kerr black hole 
in terms of gauge- and tetrad-invariant quantities. For the spin 2 black hole case the connection between the master 
equation and the de Rham Laplacian of the Riemann tensor was first noted by Ryan • Here we introduce a master 
equation for the Kerr- Taub-NUT spacetime whose symmetries allow the separation of the equation into radial and 
angular parts, generalizing some previous results valid for the Kerr and Taub-NUT spacetimes, and use it to study 
the question of superradiant scattering modes. 

II. THE KERR- TAUB-NUT METRIC 

The metric of the Kerr- Taub-NUT spacetime j3|l in Boyer-Lindquist-like coordinates (a;° = t,x^ — r, = 9,x^ = 4>) 

is 

ds^ = ^{A-a^ sin^ 0)dt^ -^[AA^a{j: + aA) sin^ e]dtd(j) 

--[{J: + aAf sin^ 9 - A^A]d(j)^ - ^dr^ - T.dO'^ . (1) 

Here E, A and A are defined by 



i:^r^ + {i + acoseY, A^r^ ~2Mr~ t + a\ A = asm^ 9 - 2icose . (2) 

Units are chosen such that G = c = 1, so that {M,a,£) all have the same dimension of length. The source of 
the gravitational field has mass M, angular momentum J = Ma (i.e., gravitomagnetic dipole moment) along the 
^-direction, and gravitomagnetic monopole moment —i. 

The two solutions r± — M ± V — a^ + P of the equation A = define the radii of the inner (r_ ) and outer 
(r+) horizons when a^ < A'P + P . Our attention will be confined to the region outside the outer horizon: r > r+. 
Adopting the Misner periodicity condition discussed by Miller , the time coordinate is assumed to be periodic with 
period Stt^, while the angular coordinates have their usual ranges: O<0<7r, O<0< 2-k. 

The metric is of Petrov type D and a Kinnersley-like null frame |19| 

I = ^[{Y, + aA)dt + Adr + ad<p] , 
n = ^[{i: + aA)dt - Adr + ad^] , 

m — —j= ^ {Acscddt — ide + c&cQd^ (3) 

V2(£ — ir + acos0) 

can be introduced to define Newman-Penrose (NP) quantities. The only nonvanishing Weyl scalar is 

^2 - {M - il)p'' (4) 

and the only nonvanishing spin coefficients are 

p = -{r - ie ~ ia cos 9)-'^ , (3 = ~ p* cot 9 / (2^2) , 

PL = p^p* A/2, TT = sin 6*7^2 , 

a — ~ (3* , T — —iapp* sin 0/\/2 , 

7 = p + pp*{r-M)/2 . (5) 
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1/2 -1/2 1 -1 3/2 -3/2 2 -2 



TABLE L The spin-weight s and the physical field component 'I' for the master equation. 



A master equation for the gauge- and tetrad-invariant first-order massless perturbations of any spin in this back- 
ground can be given starting from the following Newman-Penrose relations for any vacuum type D geometry (here 
considered with no backreaction) 



{[D-p* + e*+e- 2s{p + e)](A + p- 2s-f) 

-[5 + 7T* -a* +P- 2s(t -f P)] {S* +77- 2sa) ~ 2(s - l)(s - 1/2)^-2}* = 



for spin weights s = 1/2, 1, 2 and 



{ [A - 7* + ^* - 7 - 2s(7 + ^l)]{D-p- 2se) 

^[5* -T* +P* -a- 2s{a -|- 7r)](^ - r - 2s/?) - 2(s -)- l)(s -I- 1/2)^^2}* = 



(6) 



(7) 



for s — —1/2,— 1,-2. The case s = ±3/2 can be derived instead by following the work of Giiven ^2(?|, which is 
expressed in the alternative Geroch-Held-Penrose formalism Finally the case s = is given by 



[DA + AZ) - 8*8 - 88* + (-7 - 7* + A* + /i*)^ + {e + e* -p* - p)A 
+ {-(3* - TT + a + T*)8 + (-TT* + T - /? + a*)(5*]* = . 



(8) 



Note that only in these NP equations has the standard notation for the directional derivatives D — l^d^, A = 

and 8 = m^df^ been used and the second of these should not be confused with the equally standard notation for the 

metric quantity A — — 2Mr — + used everywhere else in this article. 

As in the case of the Kerr spacetime @, representing a gravitomagnetic dipole, and the Taub-NUT spacetime 
0, representing a gravitomagnetic monopole, all these equations for distinct spin weights can be cast into a single 
compact form in the Kerr- Taub-NUT spacetime as well, by introducing a "connection vector" with components 



1 



[M + r)a^ + {M - 3r)e^ +r^{r- 3M) 2£ + a sin^ cos 61 



p0 — _ 



- (r - M) 



air — M) cos 9 
- + i 



satisfying 



sin^ e 



1 



F^F„ 



cot^ 
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The resulting master equation has the form 



[(V^ -I- sF'')(V^ + sT^) - 4s2^2]* = 



0,±i ±1,±§,±2 , 



(9) 



(10) 



(11) 



where 1^2 is the Kerr-Taub-NUT background Weyl scalar given by Q). This master equation characterizes the common 
behavior of all these massless fields in this background differing only in the value of the spin- weight parameter s. In 
fact, the first term on its left-hand side represents (formally) a d'Alembertian, corrected by taking the spin-weight 
into account, and the second term is a (Weyl) curvature term also linked to the spin- weight value. Table |J shows the 
various Newman-Penrose quantities for which the master equation holds following the standard notation [Ij , where 
in the spin-2 case -00 and ^^4 refer to the perturbed Weyl scalars. 

The standard notation of the Newman-Penrose formalism that is employed here will be confined to this section. 
Some of the symbols of this formalism will be used to designate different physical quantities in what follows. For 
example, m and 7, used above in equations ||2J) and (jSJ, will denote respectively an angular momentum parameter 
and the Lorentz gamma factor in the rest of this paper. 
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III. SOLUTION OF THE MASTER EQUATION 

Remarkably the master equation (jllfl admits separable solutions of the form 

7Pit,r,e,(t)) = e-^"*e""'^i?(r)y(6') , 



(12) 



where w > is the wave frequency and m is the azimuthal separation constant (and eigenvalue of the usual angular 
momentum operator associated with the axial symmetry). Since the (f) coordinate has period 2tt, m must be an 
integer in order for these solutions to be smooth. Assuming the Misner condition that the t coordinate be periodic 
with period 8tt£, one is similarly led to the "quantization condition" that iuji also be an integer |^. 
The radial equation is then 



dr 



.s+idi?(r) 



dr 



(13) 



with 



V(rad)(0 = s{.s + l) + n + iu^{r^ + 2Mr + 7e^) + 



[am - 2uj{Mr + e^)f 



+2is[uj{r - M) + 



am{r - M) + 2uj{Ma^ - r{M^ + ^^)) , 



(14) 



where f2 is the separation constant. Clearly the solution R{r) of this equation depends on the value of the spin weight 
s, so when convenient this dependence will be made explicit using the notation R(r) = Rs{r). The radial potential 
can also be given a more compact form 



^(rad)(»^) 



2isuj\/ A 



dP 

dr ' 



(15) 



fl + a^u^ - s{s + 1) - 4a;2£2 - 2wam and F [r^ + + r - am/wJ/VA. 



where A 

Equation 1)13(1 will be studied on the interval r S (r+,oo), where the metric and the chosen tetrad JSJ are well 
behaved, closely following the usual treatment of black hole perturbations that motivates the present investigation. 
Of course, our analysis can be extended to the whole spacetime, but this would require introducing new coordinates 
and tetrads and would further complicate the discussion. 

By introducing the scaling 



Rir) = {r'^ +a'^ +e^)-iA-iH{r) = Qj^ H{r) 
and the "tortoise" coordinate transformation r — > , where 

A 



dr 
dr. 



+ £2 ' 



the radial equation can be transformed into the one-dimensional Schrodinger-like equation 





dr. 



—H{r)+VH{r) 



(16) 



(17) 



(18) 



with the potential 



V 



where 



(r2 + a2+£2)2 
A 

'(r2 + a2 +^2)2 
s(r - M) 



dG 2iKs{r - M) 
dr. 



(r2 + a2 + ef 



[— A + 4iwrs] 



rA _ d 1 n 

(r2 + a2 + P) (r2 +a2 + £2)2 = ^^T 2« 



(19) 



(20) 
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and K = {r'^ + + £'^)uj — am have been introduced in analogy with Teukolsky's treatment of the perturbations of 
the exterior Kerr spacetime. It is useful to rewrite this potential in the following more compact form 

y ^ _n-i^Q [A-^-2»X-s(r-M) + A(-A + Wg)] 

The asymptotic form of the radial equation as r — > oo (r.^ ^ oo) is 

^H{r) + (c.2 + ^) H{r) ^ , (22) 

which has asymptotic solutions H ^ j,±Sg^jwr»^ ^ ^ ^-lujr, ^^^^ ^ ^lujr^ j^is+i ^-^^ accordance with the 
peeling theorem and in this regime the effect oil appears negligible. On the other hand close to the horizon 

r ^ r+ (r, ^ — oo), the asymptotic form of the radial equation becomes 



^i/(r) + (fc-ife+)'i/(r) = 0, (23) 



where 6+ = 2(Mr++^^) i^ — ^^ r7ia;+ and 

o a . 

" r2 + a2 + £2 = 2(Mr+ +£2) ( ) 

is the "effective angular velocity" of the horizon. The asymptotic solutions are H ^ g±»(fc-«f)+)''. ^ ^±s/2g±jfcr, ^ j 
]l ^ e*'^'"' and R ^ ^-ikr, ^ Only one of these two behaviors for R is correct in the sense that it implies regularity 
of the fields on the horizon r = r_|_ as stated by Teukolsky. The asymptotic boundary conditions on the horizon are 
a delicate problem because both the coordinates and the tetrad being used in this discussion are singular there, but 
it was solved by Teukolsky in the Kerr case by picking out the second solution as having the correct behavior in that 
limit, which must therefore also be the correct choice in the present more general case. This result follows from the 
requirement of causality that at the horizon we must choose the boundary condition that the wave is always ingoing. 

An immediate consequence of the behavior of the solution at the horizon and at spatial infinity is that there are 
superradiant scattered modes as in the Kerr case jlll l22j , but now infiuenced by the nonzero value of the parameter I. 
For superradiant modes, reflected waves carry away more energy than the incident waves bring in, with the rotation 
of the black hole supplying the extra energy. 

To understand how this comes about, consider first the 5 = case and imagine a solution of equation l|18|) cor- 
responding to the reflection and transmission of a radially incident wave of frequency a; > 0. For — )■ oo one has 

$~^e-*"(*+''*)+y^7ee-''^(*-''*\ (25) 
where A is the amplitude of the incident wave and is the reflection amplitude, while for 7'* ^ — oo one has 

$-yire-*"*-*'=''% (26) 

where T is the transmission amplitude, and the angular coordinates in $ are suppressed. The potential Y is real in 
this case (s = 0); therefore, it follows from flux conservation that 

|7^|2 + -|r|' = 1. (27) 

LlI 

If fc/w < 0, then |7?.| > 1 and so one has superradiance. The superradiance condition 

k = bj- mcj+ < (28) 

thus depends on the value ujj^ of the effective angular velocity of the horizon. For any frequency w > 0, this condition 
can be satisfied for large enough values of the azimuthal separation constant m of the same sign as w+ and a, i.e., 
when the angular momentum of the wave is in the same sense as the angular velocity of the horizon. It follows from 
Eq. (|24|l that, for fixed a = J/iV/, r^. increases and hence lOj^ decreases if M or \l\ are increased. Thus for a spherical 
system where a = implies lo^ = 0, there is no superradiance, and in the KTN case where a 7^ 0, this phenomenon 
can also be suppressed by increasing M and/or \l\ such that u;+ becomes so small that, for fixed u and m, the quantity 
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k becomes positive. In fact superradiance can be suppressed independently of the values of M and J by making \£\ 
sufficiently large. Note that the transmitted wave in Eq. (|26l) moves toward the horizon for k/oj > 0, but reverses 
direction in the case of superradiance with k/oj < 0; moreover, total reflection occurs in the exceptional case that 
k/uj = 0. 

When the spin-weight s is an integer, we would expect that there would still be superradiance because in the ^ — > 
limit one obtains the Kerr case which is well known to exhibit this phenomenon for integer-spin fields. In the case 
of half-integer spin fields, the Kerr case does not exhibit superradiance, but when £ ^ this effect could in principle 
change, so the question must be re-examined which we do here using a field theoretic approach 14, 20, 23]. In fact one 
condition for having superradiance is that the flux of particles across the (null) horizon in the forward time direction 
be negative 

{^^Ji:^)\r=r^<0 , (29) 

where S. = dt + oj+d^ is the future-directed (null) normal to the horizon and J|'^| is the conserved particle number 
current vector associated with the various fields of spin \s\. This condition corresponds to a positive flux out of the 
horizon. 

Alternatively one may consider the rate {dN/dt)i-a at which particles are falling in through the horizon per unit 
time, which must be negative for superradiance to occur 

where g is the determinant of the spacetime metric. The appendix shows how this integral relation follows from 
integrating the particle flux over a suitable region of the spacetime. 
For the scalar field case one has 

Jo^ =in($*V^$-$V^$*) , (31) 

so that 

{i,Ji;)\r=r^=2hk\A\''\n\Ym\ (32) 

which is negative when k — lo ~ muj^ < 0, corresponding to superradiance. Analogously, evaluating the number of 
particles entering the horizon per unit time gives 

^) =2hk^\Am\ (33) 

dt J i„ UJ+ 

which is negative under the same condition. Note that the angular part Y{9) of the scalar field, which can be taken 
to be real, has been normalized by 



/ Y^{0) siTi0d9 = — 
Jo 27r 



(34) 



In fact since A = on the horizon, the surface element for a sphere with the horizon radius is just + aA)\r=r+'i9d(j), 
where (E -I- aA)\r=r+ = r"^ + + l"^ ^ so the integral of (|32|l over the angular variables produces exactly (|33|l . 

In contrast with the case of bosons, fermions do not exhibit superradiance, mirroring exactly the corresponding 
behavior in the Kerr spacetime. For the \s\ = 1/2 case, the neutrino particle number current is 

^ ^^P^ ^ (35) 

where 7'' are the (coordinate) Dirac matrices. Passing to the more standard spinor formalism (see |l4| p. 539 for 
notation and conventions), this current becomes 



V2' 



(36) 



with 



-P^ and g 



A' 



P^' and 



V2 



If" mf" 



(37) 
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In particular one finds 

= -XI , = Xo , = -Xl , P'' = Xo ■ (38) 

Since the spinor formalism can be immediately converted into the Newman-Penrose language, the neutrino current 
can be also be expressed as 

^1/2^ = 2r\Xi\' + n^\Xo\' - m^(xixS) " ™*^(XoxD] • (39) 

In this case (|s| = 1/2) the condition for superradiance is never satisfied. In fact following Giiven ,20,] and using the 
relations 

IXol' = |i?i/2(0Pn%: IXiP = ^\R-i/2ir)\'Y\/,, (40) 

in (jSnj and (ISOJ, one finds 

- (A|i?i/2(r)P - 2|i?_i/2(r)n 1.=.+ . (41) 
Next using our previous results for the behavior of Rsir) near the horizon 

\Ri/2{r)\' - \R-i/2{r)\' - |c_i/2pA (42) 



with c-i-i/2 constants, one finds the final result 



|ci/2p>0 (43) 



which means that there is no superradiance in this case. 

In the Rarita-Schwinger case (|s| = 3/2), the field is described by the Majorana spinor-valued 1-form 0113 
satisfying the further conditions V^^" = and 7q^" = 0, and the current is 

J3/2'' = ^al'^'i" . (44) 
We expect that a calculation similar to that valid for the neutrino field would show that 

e^'I'a7''*"|r=r+ < , (45) 

namely that even in this case there is no superradiance. 

Of course the study of the equations for the remaining (gauge- and tetrad-dependent) components of the various 
fields (i.e. 0i, ips, etc.) and of theorems analogous to those found by Wald [23|, Fackerell and Ipser and Giiven 
|20l | for the Kerr case require some extensions here. For the Kerr spacetime these theorems essentially state that the 
master equation is enough to describe all of the relevant physics of the fields considered in the exterior spacetime. 

In the Kerr-Taub-NUT case the master equation completely describes the scalar (s = 0) and neutrino (s = ±1/2) 
fields over all of the spacetime outside the horizon. For higher spin fields of course the master equation describes only 
the highest and lowest spin- weighted components, as specified in Table 1. However, because of the "peeling theorem" 
|l3l I23 , this equation is sufficient to understand the relevant physics of the fields of any spin considered at spatial 
infinity, in the sense that other spin-weighted components become negligible. Close to the horizon, one has instead 
only a subset of perturbations under control, namely those indicated in Table 1. 

To complete our analysis of the perturbations on the Kerr-Taub-NUT spacetime, we need to discuss the angular 
equation 



1 d 



sin6'd6l V d6' y 

where 

(s + 2^w)2m cos + {s + 2£lu)'^ + 



+a^Lu'^ cos^ + 2auj{-s + 2euj) COS 9 . (47) 
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This equation generalizes the spin-weighted spheroidal harmonics of Teukolsky^J |28L l29l |. Here instead of the 
parameter s in the spin-weighted spheroidal harmonics, the combinations s + 2£uj and —s + 2£ll! appear. In fact, 
introducing a new variable x — 1 + cos 9 and the rescaling 



one gets 



where 



Y = a;l™-'^-2c^'^l/2(2 _ a;)l'"+''+2'^^l/2x(a;) , (48) 
x{x - 2)X" + {Bix + Bq)X' + {C2X^ + Cxx + Cq)X = , (49) 



Bo = -2(|TO-s-2tj^| + l) , 

Bx = 2 -f |to - s - 2uji\ + \m + s + 2uj£\ , 

Co = ^ [m^ + {s + 2ujlf + - (s -\- 2^if\ + |m + (s + 2ijjl)\ + |m - (s + 2uj()\\ , 

+r2 - aLj[acj + 2(s - 2cj^)] , 
C\ = 2acj(aa; -|- s — 2a;i?) , 

C2 = -a^w^ ^ (50) 
which is a generalized spheroidal wave equation of the Leaver form . 

IV. GRAVITOMAGNETISM 

We note that these results contain the combinations ±s -f 2ujI^ so that, in a certain sense, the spin is coupled to the 
gravitomagnetic monopole moment; for a = only the combination s -\- 2lo(. is involved. This is a novel manifestation 
of the spin-gravity coupling jSl], . 

To interpret this coupling let us start with the propagation of test electromagnetic fields in the KTN spacetime. 
Using the Skrotskii formalism |3.l| | in which the gravitational field may be replaced by an equivalent material medium, 
the problem reduces to Maxwell's equations in a certain gyrotropic medium in a global inertial frame with Cartesian 
coordinates (t, a;, y, z). In this section only we employ the other metric signature — h which is more standard for 
the following considerations. To simplify matters we assume that the electromagnetic waves have a time dependence 
of the form e~''^*; moreover, we linearize the KTN metric. With these simplifications, the wave equation becomes 



where P = -zfiV, W=t = E ± iH, 



[P + 2fta;A(g)] x W± =F«ftw7VW± , (51) 



3 X p z, X p , , 

Afe) = -TTT + , 9 , 52 
is; |p|3 ^2 ' V y 



and A/" ~ 1 + is the index of refraction of the medium. Here E and H arc complex fields and \p\ — y/ x"^ + y'^ + z'^ 

is the isotropic radial coordinate in the linearized KTN spacetime. Furthermore, W"^ represents the positive helicity 
wave amplitude and W~ represents the negative helicity wave amplitude. It is a general result that helicity is 
conserved in pure gravitational scattering. A more detailed treatment is available elsewhere [s^l ■ 

As stated above, a satisfactory interpretation of the KTN spacetime requires that the time coordinate t be periodic 
with period 8tt£. Thus a propagating wave with a time dependence of the form e~*'^* is possible only when Aui is 
an integer. Eq. (|51|l can be expressed in the Dirac form and then interpreted as the wave equation for a particle of 
inertial mass mo — huj and gravitomagnetic charge qo = —2fiw propagating in a gravitomagnetic field B(g) = V x A(g) 
given by 

B(g) = |^[3(p-J)p-J]-^p^. (53) 

Restricting our attention now to the motion in a purely monopole field (i.e. J = 0) of strength /io = — ^, wc note 
that the classical equation of motion would be 

p X V 
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where p = mo7V is the kinetic momentum. It is well known that this equation has a constant of the motion given 
by = p X p — goMoPi which can be interpreted as implying that the orbital angular momentum of the particle is 
augmented, through its interaction with the monopole, by S' = —qonop = 2hu!£p. Note that the magnitude of this 
vector S" = |S'| is a positive integer multiple of fi/2, since Acol is an integer. The intrinsic rotational symmetry of 
the monopole thus leads to a conservation law for the total angular momentum of the particle that consists of 
an orbital part plus a contribution from the angular momentum of the total field (generated by the monopole and 
the particle) that is acquired by the particle through its interaction with the monopole. That this mechanical result 
extends to the wave treatment of the scattering problem has been demonstrated by a number of authors in the case 
of magnetic monopoles (see e.g. [s^H^). 

Note that the canonical momentum in the present case would be P = p + A(g) = p — 2/ia;A(-g) . The total angular 
momentum for a spinning quantum particle turns out to be 

J = px[P + 2hujAi^g^] + S + S' = L + L' + S + S', (55) 

where L = pxP (satisfying [Li, Lj] = iheijk L^) is the orbital part, S is the spin part, L' = 2hujpxAi^g^ — —2huj£ cot 90 
and (satisfying [J7i, JT}] = iheijk <Jk) is the generator of spatial rotations. Hcuristically one may say that the spinning 
particle, in its interaction with the monopole, picks up an additional spin contribution S' — 2hwlp, such that the net 
effective spin of the particle along any radial direction (like the z— axis) — taken to be the axis of quantization — would 
be S + S' . . , —S + S' . The general situation is, however, not so straightforward since ^ = L + L' + S + S' and we 
note that 

[U + L[,Lj + L'j\ = ihe,jk{Lk + L'k + S'^) . (56) 

Thus the separation of the total angular momentum H55I) into orbital and spin parts is not quite obvious in this case. 
Nevertheless, since the master equation describes the behavior of the highest and the lowest spin-weighted amplitudes, 
it is natural to expect that the angular part of the equation would only contain the spin-weight combinations s + 2uj£ 
and — s + 2uji. 

V. DISCUSSION 

A master equation for the gauge- and tetrad-invariant first-order massless perturbations of any spin s < 2 on the 
Kerr-Taub-NUT background spacetime has been obtained and separated. We have studied superradiance in this 
case and have shown that the situation is very similar to the Kerr spacetime; in particular, we have demonstrated 
the absence of superradiance for half-integer spin perturbations. Furthermore, the interaction of the perturbing field 
with the gravitomagnetic monopole contributes a certain half-integer spin component to the angular momentum 
that combines with the spin of the field; this novel form of spin-gravity coupling has been briefly discussed here. 
This investigation offers the possibility of achieving a better understanding of perturbations of black hole spacetimes 
within this larger family. That is, one can in principle extend the methods developed in [3^ to discuss various 
aspects of scattering of radiation from the Kerr-Taub-NUT spacetime such as the polarization properties of the 
scattered radiation, glory effects and quasinormal mode oscillations. These may then be used to search for rotating 
gravitomagnetic monopoles. 

The source of the Taub-NUT solution is a gravitational dyon (with gravitoelectromagnetic monopoles m and 
—£). The Taub-NUT solution can be extended to include a cosmological constant or an infinite set of multipole 
moments pertaining to axisymmetric deformations of a rotating source . The periodicity of the time coordinate 
in the Taub-NUT spacetime renders such sources, if they exist at all, rather exotic astrophysical systems. Normal 
astronomical systems do not exhibit such periodicity; therefore, the temporal periodicity may be associated with 
Hubble-scale structures [s^ or compact dark-matter candidates in these connections, lensing properties of Taub- 
NUT spacetime have been studied in detail [3^ ^3 and compared in the latter case with astronomical observations 

A remark is in order here concerning an indirect application of Kerr-Taub-NUT spacetime to rotating relativistic 
disks, which are of great interest in astrophysics. The KTN spacetime can be employed to generate exact solutions of 
the Einstein-Maxwell equations corresponding to stationary axially symmetric disklike configurations of matter with 
a magnetic field ,41^, . 

The Taub-NUT spacetime has had significant applications in theoretical studies of the spacetime structure in 
general relativity and, more recently, in quantum gravity. Its Euclidean extension is important for the study of 
monopoles in gauge theories. Embedding the Taub-NUT gravitational instanton into five-dimensional Kaluza-Klein 
theory leads to a Kaluza-Klein monopole (ijj. [il . . Euclidean Taub-NUT spaces have been discussed by a number of 
authors in connection with monopoles in supersymmetric gauge theories 46] . Further generalizations and extensions 
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of Taub-NUT spaces (such as the Kerr-Newman-Taub-NUT-adS spaces) are topics of current research in string theory 



APPENDIX A: INTEGRAL FORM OF THE PARTICLE NUMBER CURRENT CONSERVATION LAW 

This article has dealt exclusively with the exterior KTN spacetime. Therefore, in the determination of the number 
of particles falling in through its horizon per unit time, one needs to apply the conservation law for the particle number 
current only in the exterior region. To this end, consider the volume enclosed by two "spheres" of time-independent 
radii r-|_ < i?i < i?2- Eq. H30|l can be derived easily from the integral of the particle number conservation equation 

= V^J";a = 9a (\/^ J") = dt{V^j') + 8,(^1') (Al) 

over the region V within a time coordinate hypersurface bounded by these two spheres, and then using Gauss's law 
to convert the second term to a surface integral over the boundary 



dNv _ d 
dt ~ di 



[ ^J'drded(t) = - [ d,{y/^.r)drded(l3 
Jv Jv 



^f^rdOdct) + / y/^,rd9d<j) . (A2) 

r=n.2 Jr=Ri 

This gives the rate of change of the number Ny of particles in this region in terms of the flux entering the outer 
sphere and exiting the inner sphere. 

When the integration domain is replaced by the region enclosed between the horizon and a large sphere at infinity, 
one may compute the rate at which particles are leaving this region through the horizon alone. The contribution from 
the horizon, i.e., 

'^rdedci) = ^ + / y/^rded<j) , (A3) 



r— r+ ^ r—oc 

when sign-reversed, gives the rate at which particles are entering the horizon, leading to Eq. I|30|) . This sign-reversed 
quantity must be negative for superradiance to occur. 
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